Impurity-induced magnetic order in low dimensional spin gapped materials 
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We have studied the effect of non-magnetic Zn impurities in the coupled spin-ladder 
Bi(Cui_xZnx)2P06 using ^^P NMR, muon spin resonance (/iSR) and Quantum Monte Carlo simu- 
lations. Our results show that the impurities induce in their vicinity antiferromagnetic polarizations, 
extending over a few unit cells. At low temperature, these extended moments freeze in a process 
which is found universal among various other spin-gapped compounds: isolated ladders, Haldane or 
Spin-Peierls chains. This allows us to propose a simple common framework to explain the generic 
low-temperature impurity induced freezings observed in low dimensional spin-gapped materials. 



The physics of low-dimension quantum antiferromag- 
nets (AF) is fascinating and surprising. Simple AF spin- 
chains or ladders display exotic low temperature be- 
havior such as spin-liquids, spin-gaps, magnetic orders, 
spin glasses, etc. The precise behavior depends on the 
value of the spin, the dimensionality of the material, the 
anisotropy, the relative strengths and signs (possibly frus- 
trating) of the magnetic couplings. Studying the effect 
of impurities is an efficient way to classify and reveal the 
quantum nature of these systems [1]. For example, the 
qualitatively different properties of half-integer and in- 
teger spin chains are evidenced in the strongly different 
impurity-induced effects in the two systems. 

In AF chains and ladders, non-magnetic impurities in- 
duce paramagnetic clouds made of alternating moments 
in their vicinity, the shape and extension of which is 
directly linked to the type of electronic correlations in- 
volved in the ground state. At low enough temperatures, 
impurities may even lead to magnetic order or spin freez- 
ing, a sort of "order by disorder" phenomenon [2, 3, 4, 
5, 6, 7, 8]. For example, a few % of non-magnetic Zn or 
Mg substituted at the Cu site of the spin-Peierls chain 
CuGeOa [3, 7, 8], the spin-ladder SrCu203 [4, 9] or at 
the Ni site of the Haldane chain PbNi2V208 [5, 10], leads 
to a collective freezing below a 3D ordering temperature 
Tg of a few K. The induced paramagnetic clouds can 
be viewed as effective localized moments extending over 
a finite region of typical size ^ around each impurity, 
with an exponentially decaying staggered magnetization 
{S%r)) ~ (-!)'• exp(-r/0 [11, 12, 13]. A 3D effective 
interaction between them governs the freezing at low tem- 
perature. The frozen-state characteristics should then be 
specific to a given compound and geometry, explaining 
why no common picture has been given yet to under- 
stand the numerous experiments. As the effective inter- 
action between impurities falls off exponentially with the 
distance r [14, 15, 16, 17] as J°^{r) - Joexp(-r/^), it 
has then been suggested that the relevant energy scale 
for the magnetic freezing is given by the typical coupling, 
i.e. the effective interaction occurring at the average dis- 
tance between impurities along the chains or the ladders 
(r) ~ 1/a; [7, 15, 18] where x is the impurity concentra- 



tion. However, as pointed out in Ref. 19, realistic pa- 
rameters would give exponentially small Tg, and a quan- 
titative agreement with experimental values has been ob- 
tained using unrealistic enhanced lengths ^ [7, 18]. 

To decide which actual process drives the freezing in 
impurity-doped spin-gaped materials, we have chosen to 
study the effect of impurities in a new spin ladder mate- 
rial BiCu2P06 (BCPO). In contrast with the archetypal 
ladder compound SrCu203 where ladders are almost iso- 
lated with a magnetic coupling Jiog ~ 2000 K, BCPO has 
a non- negligible inter-ladder coupling but a much smaller 
Jicg 100 K [20]. We present here ^SR and NMR stud- 
ies of the effects of non-magnetic Zn impurities, together 
with Quantum Monte Carlo (QMC) simulations. Above 
Tg, NMR enables one to measure the induced param- 
agnetic cloud near the impurities while /iSR is an ideal 
probe of the frozen state, measuring both the freezing 
temperature Tg and the corresponding field distribution 
in the full volume of the sample. Our results compared 
to QMC simulations demonstrate that the extension of 
the induced clouds is only a few unit cells. This small 
extension, together with the 3D coupling J3D, are the 
crucial parameters in the freezing process. We thus ex- 
plain in a common framework the freezing temperatures 
observed in ladders and other low dimension spin-gapped 
materials. 

Synthesis and structural characterization of 
Bi(Cui_xZnx)2P06 are presented elsewhere [20, 21]. 
NMR measurements were performed using a 7 Tesla 
homemade spectrometer with standard pulse techniques 
and Fourier Transform recombinations. Magnetic 
susceptibility was measured in a SQUID magnetometer 
at H=0.1 T. When Zn is substituted at Cu site, it is 
expected to release a free 5* = 1/2 spin. Macroscopic 
susceptibility indeed displays the signature of this 
free spin through a low temperature Curie component 
proportional to Zn content. The fitted Curie constant 
corresponds to a total induced spin S = 0.35 — 0.45 per 
Zn slightly smaller than the expected = 1/2 because 
of possible frustration effects [21]. This total spin is not 
localized on a single site but develops as an induced, 
extended, alternating cloud (AF cloud) along the ladder 
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FIG. 2: (Color online). QMC results for the AF cloud ex- 
tensions and the correlation length of model (1). (a) T de- 
pendence of the induced extensions and for coupled and 
isolated ladders (SrCu203 has Jicg/Jrung = 2). (b) T = 
results for the inverse gap J/ A ~ ^ and the extensions S^x, 
versus J±. The black line is the pure correlation length 
^ ~ (Jx — Ji_)~'^ (see text). The small error bars come from 
the uncertainty in the 3 parameter fit of QMC data to Eq.(2). 



FIG. 1: (Color online). Lower panel: Left: ^^P NMR spectra 
for BCPO: Zn 2% at H=7 T; Right: QMC results for the 
coupled ladders model Eq. (1) at T = J/10 with an applied 
field H = J/10, using Jiog = Jrung = J, J±/J ~ 0.1, x = 
2%. The corresponding Zn induced magnetic pattern in the 
ladders is displayed in the top panel, with a S=l/2 arrow 
indicating the scale. 



and in adjacent ones. This AF cloud reveals itself in 
the ■^^P NMR spectrum through a low temperature 
broadening (Fig. 1) just as in isolated ladders [9]. The 
spin-i coupled ladders can be modeled by: 

+ J±Si^2j + l ■ Si^2j+2- (1) 

As already studied in details [24, 25], this model dis- 
plays a gapped valence bond solid (VBS) phase when 
the inter-ladder coupling J± is not too strong. For 
isotropic ladders (Jrung = ^icg) this is the case when 
Jj_ < Jl with Jl/J = 0.31407(5) [24], and the corre- 
lation length diverges close to the quantum critical point 
(QCP) C {Jl-J±y with u = 0.709(6) [25]. Previous 
bulk studies on BCPO [20] estimate Jiog ~ Jrung '^-^ 100 
K. Muffin tin orbital calculations suggest the existence of 
an additional frustrating coupling along the ladders [20] . 
We did not take this possible frustration into account 
in our model Eq. (1) as we just want to capture semi- 
quantitative features of the low-T physics. The experi- 
mental spin gap A ~ 35K yields an inter-ladder coupling 
0.1 < Jj_/Jieg < 0.2, implying a small C Using QMC 
simulations we have computed the magnetization profile 
induced by a single impurity on a 32 x 32 lattice. Ex- 
ponentially localized AF profiles are found, with the 2D 



form 

(5^(r)):^5oexp(-f -f)(~l)-+'' (2) 

valid all over the gapped regime. The extensions and 
^j, at r = are reported in Fig. 2(b) where one can 
clearly see that they quantitatively follow the behavior 
of the pure correlation length ^ ^ J/ A: the impurity in- 
duced effect reveals the intrinsic properties of the ladder, 
as in Haldane chains [22, 23] or high Tc cuprates [26]. 
This result contradicts analysis of Ref. 18 which argues 
that the AF cloud extension depends on the impurity 
content and reaches values as large as 50 cell units at low 
concentrations. However, such conclusions were entirely 
based on NMR experiments done at the very low x= 0.1% 
concentration where the observed NMR broadenings are 
too small to be reliable. All other data of Refs. 18 and 
9 are fully compatible with our own results and conclu- 
sions. As shown in Fig. 2(a), the cloud extensions and 

increase when the temperature is lowered but converge 
very rapidly once the gap A ~ 0.2 — 0.3 J is reached. As a 
result, the NMR broadening (Fig. 1) at low T is expected 
to follow the Curie law due to 5*0 only. This is indeed the 
case in the experiments presented in Fig. 1. QMC sim- 
ulations performed for J±/ J = 0.1 at T = J/ 10 with an 
external field H ^ J/IO over 32 x 32 sites with x= 2% of 
non-magnetic impurities are averaged over 10^ disordered 
samples to obtain the distribution of local fields seen by 
each phosphorus plotted in the right panel of Fig. 1. The 
simulated broadening found without any adjustable pa- 
rameter is remarkably similar to the experimental ones, 
albeit possible frustration is not considered. 

In order to probe whether these impurity induced AF 
clouds of only a few unit cell extension lead to a freez- 
ing in BCPO, we performed a /iSR study on the same 
Zn substituted BCPO samples at the PSI (GPS) facil- 
ity. Pressed disks of randomly oriented powders were 
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FIG. 3: (Color online). Muon polarization for BCPO above 
and below Tg (black circles) on upper and left axis is com- 
pared to a similar experiment made in SrCu2 03 (bottom axis 
being rescaled, right axis is slightly shifted because of a dif- 
ferent background) from Ref. 27. Grey line is a fit to Eq.3 of 
Ref. 27 as explained in the text. 
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FIG. 4: (Color online) Left panel: transition temperatures 
versus impurity concentration for various low-D spin-gapped 
systems: coupled ladders Bi(Cui_a;(Zn or Ni)a;)2P06 from 
this study; isolated ladder Sr(Cui_a;(Zn or Ni)^)203 [4] [18]; 
Haldane chain Pb(Nii_i,Mgi:)2V208 [31]; spin-Peierls chains 
Cui_i(Zn or Ni)a;Ge03 [8]. Full and open symbols corre- 
spond respectively to non magnetic and magnetic impurities. 
Right panel: Same data where Tg is rescaled by its value at 
x=3%. 



used. Muons which carry a spin 1/2 are implanted in 
the samples and precess around the local magnetic field. 
This precession is measured through the detection of the 
asymmetry of the positron emission due to the muon de- 
cay. In the absence of an external magnetic field ("zero 
field setup"), above Tg, only the tiny nuclear spin dipo- 
lar fields are experienced by the muons. This results in 
a gaussian field distribution and a very slowly-decaying 
" Kubo-Toyabe" polarization as displayed in Fig. 3 at 
T=10K. When the electronic spins start to freeze be- 
low Tg, their randomly oriented static moments result 
in a much larger field distribution and a much faster de- 
caying asymmetry as shown in Fig. 3 at T=1.6 K. We 
checked that this depolarization is of static origin below 
Tg by applying a large longitudinal field and finding the 
expected asymmetry decoupling. Since no marked oscil- 
lations are observed, the corresponding field distribution 
is not that of a simple commensurate magnetic ordering. 
However, due to the lack of a dip in the asymmetry, it is 
probably not completely random as in a spin glass, but 
more likely in an intermediate situation with AF clusters, 
as discussed in Ref. 27. To get an estimate of the field 
distribution, we used the same phenomenological expo- 
nential field distribution as in Ref. 27, which fits well the 
data as shown in Fig. 3. The resulting static field distri- 
bution develops as a mere order parameter which allows 
to determine Tg. Whatever the impurity nature or con- 
centration, the transition temperature is proportional to 
the field distribution even for very large impurity con- 
centrations. 

The SrCu203 ladder shows similar time dependence 
for the muon polarization as displayed in Fig. 3 [27]. 



Both the static field distribution and Tg [4] are found 1.7 
times larger than in BCPO. This 1.7 ratio is unexpect- 
edly about one order of magnitude smaller than the ratio 
between the ladder couplings. Other low-D spin gap sys- 
tems such as Haldane chains or Spin-Peierls chains also 
display very similar Tg despite their different geometries 
and gap values (Fig. 4). The dependence on impurity 
content also shows a strikingly generic behavior, with a 
similar departure from a linear behaviour at about 2-3% 
for all compounds as demonstrated in the right panel of 
Fig. 4. Magnetic impurities in these materials also lead 
to similar Tg and similar concentration dependence as 
shown in Fig. 4. 

Despite several theoretical investigations devoted to 
understand the origin of the impurity-induced 3D or- 
dering in low-D gapped systems [15, 17, 19, 29, 30], no 
common framework has emerged so far which explains 
this generic Tg behavior. The collective freezing of the 
effective moments (having a 3D extension ~ S.xS.y^z at 
T > Tg) is actually controlled by the exponentially de- 
caying 3D coupling of the general form 

\JlE{f)\ ^ Jaoexp - |- _ , (3) 

expected to occur for the wide class of spin gapped ma- 
terials [14, 15, 16, 17, 32, 33]. Couplings in the three 
directions are necessary to allow finite-T ordering due to 
the Mermin- Wagner theorem, and J3D is the weakest of 
these couplings. The typical coupling Jtyp = I >^3d ( ) I ; 
taken at the average distance (r) between impurities, is 
much too small to explain the actual 3D freezing temper- 
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atures [19] because it does not take into account the rare 
but crucial situations where r is small. On the contrary, 
the average coupling Jav taken over all possible J§^{r) 
does account for the broad distribution of effective inter- 
actions and is just given by 



^av=(|J3°g(r)|)=. J; 



3D 



(4) 



where ^ ^x£,y£,z is the magnetic volume occupied 
by each induced moment. We propose that this aver- 
age coupling governs the ordering, i.e. Tg ~ Jav In- 
deed, this model accounts well for the fact that all Tg 
are similar in magnitude, because both J3D and are 
roughly of the same order of magnitude among the dif- 
ferent systems. Precise values for J3D are not known 
since they are at most a few percent of the dominant 
coupling. However, rough estimates give a few K for 
CuGeOa [34] and PbNiiViOg [5]. For SrCuaOg a theo- 
retical calculation [29] based on an isotropic ladder model 
found J3D 20K. Even though the order of magnitude is 
correct, the actual J3D is probably smaller because of the 
anisotropy Jieg/Jrung = 2 [28]. Regarding BCPO, based 
on our analysis we can confidently predict a 3D coupling 

5K. Turning now to the magnetic volumes of the 
induced moments, they are similar at low temperature, 
typically ~ 30 sites, despite the geometric differences be- 
tween systems. For instance, Zn in the almost isolated 
anisotropic ladders Sr(Cui_a;Zn2;)203 induces a cigar-like 
pattern, with ~ 2 x 8 x 2 ~ 32. In the coupled lad- 
ders Bi(Cui_j;Zna;)2P06 with J±/J ~ 0.15, Zn induces 
a pancake-like volume ~ 2 x 3.9 x 2 x 2.1 ~ 33. This 
is because all these spin-gapped materials are far enough 
from a QCP so that the various ^ remain a few unit cells 
and do not diverge (Fig. 2), whereas a totally different 
physics would be observed at criticality [35]. 

Such an analysis allows to understand the concentra- 
tion dependence of Tg as well. At small x, the average 
volume between impurities (r)'^ ^ is large compared 
to the magnetic volume V^, and Eq. (4) implies a lin- 
ear behavior of Tg ~ Jav with x. At larger x, extended 
moments get close enough to percolate and the critical 
temperature is not linear with x anymore. This should 
occur at X ~ 1/V^ ~ 1/30 ~ 3%. This prediction of a lin- 
ear increase of Tg below x ~ 3% and a deviation above 
is exactly what is observed in Fig. 4 (see also the at- 
tached supplementary material where QMC simulations 
are presented). 

In summary, the framework we propose to explain the 
impurity induced freezings in various systems is based on 
the fact that the impurity induces an extended magnetic 
moment, the volume of which is similar in all systems be- 
cause all lie far from the QCP. These induced moments 
interact through a 3D interaction which is not much sen- 
sitive to the geometry or the gap nature. This holds 
for magnetic impurities as well. It will be interesting 
to understand if this universal picture survives when a 



Quantum Critical Point is approached. 
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